Spin tunneling via dislocations in Mni2 acetate crystals 
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We show that dislocations should be the main source of spin tunneling in Mni2 crystals. Long-range 
strains caused by dislocations produce broad distribution of relaxation times that has been seen in 
many experiments. When the external magnetic field is applied along the c-axis of the crystal, local 
rotations of the magnetic anisotropy axis due to dislocations result in the effective local transverse 
magnetic field that unfreezes odd tunneling resonances. Scaling law is derived that provides universal 
description of spin tunneling for all resonances. 

PACS numbers: 75.45. +j, 75.50.Tt 



Mni2 acetate crystals exibit quantum magnetic phe- 
nomena at a macroscopic scale. They have centered 
tetragonal structure with a = 17.319A and c = 12.388A 
as lattice parameters [Gl^ . Spin-10 Mn 12 molecules at the 
sites of the lattice show magnetic bistablity due to the 
65K barrier between the spin-up and spin-down states || . 
Quantization of spin levels manifests itself in a spectacu- 
lar stepwise magnetic hysteresis 0. Other observations 
include memory effects non-exponential relaxation 
|B0 , and a peculiar crossover between thermal and quan- 
tum behavior Theoretical models explain tunnel- 
ing in Mni2 by phonons [^|Jg| Jl0| , pd]| , nuclear spins ||],|l3| , 
dipolar fields US fourth 

-order magnetic anisotropy 
||, Landau-Zener effect HHH], Jahn-Teller effect @, 
etc. Explanation of some observations remains contro- 
versial though. 

Firstly, there is no agreement among researchers on 
what causes tunneling in Mni2. The forth-order trans- 
verse anisotropy cannot account for odd resonances. Be- 
sides, it is too weak, as are hyperfine and dipolar fields, 
to provide the actual tunneling rate. Secondly, in all 
samples studied to date, a "minor species" of Mni 2 has 
emerged that exhibits faster magnetic relaxation than 
the "major species" |l9). Finally, the time dependence 
of the magnetic relaxation in Mni2 has not been under- 
stood. If all Mni2 molecules were subject to the same 
crystal field, the magnetic relaxation would be strictly 
exponential in time. In reality, however, the relaxation 
deviates from the exponential in the kelvin range and is 
clearly non-exponential in the subkelvin range EM] . 

In this Letter we show that dislocations, undoubtedly 
present in Mni2 crystals, provide long range deforma- 
tions (see Fig. 0) which should be the main source of 
spin tunneling. The "minor" and "major" species differ 
on the strength of the deformation, that is, on their dis- 
tance from the dislocation cores, with no sharp boundary 
between the two. Local rotations of the anisotropy axes 
due to dislocations are responsible for odd tunneling reso- 
nances. Broad distribution of deformations causes broad 
distribution of tunneling rates. We compute the relax- 
ation law in a crystal with dislocations and show that it 



obeys scaling that can be tested in experiment. 
We study Hamiltonian 
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where S z is the z-component of the spin operator, S = 10, 
D = 0.65 K, H z is the magnetic field applied along the 
z-axis (c-axis of the crystal) , and TL mc is the magnetoe- 
lastic coupling. The Hamiltonian of Mni2 also contains 
crystal fields of fourth order on the spin operator, mag- 
netic dipole interactions, and hyperfine interactions. We 
neglect them in order to emphasize the effect of dislo- 
cations. The magnetoelastic coupling in Mni2 is of the 
form (l|@ 

Wmc = 9iD(e xx — Ey y )(S x — Sy) + g2De xy {S x , S y } 
+ 9zD(e xz {S Xl S z ] + e yz {S y , S z }) 
+ g 4 D(aj xz {S x , S z } + uJ yz {S y , S z }), (2) 



where 
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are linear deformation tensors, u being the displacement. 

The full classification of relevant elastic deformations 
due to different types of dislocations will be done in a 
longer paper. In this Letter, we will illustrate the effect of 
static deformations on tunneling by considering edge dis- 
locations running perpendicular to the anisotropy direc- 
tion along the y-axis. For the YZ extra crystallographic 
plane inserted at z > one obtains (see, e.g., Ref. po|): 
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where < a < 1/2 is the Poisson elastic coefficient (we 
will use a = 0.25). If the extra plane is inserted at z < 0, 
the above expressions change their sign. 



Transverse anisotropy: Edge Dislocation jjY, g 1 = g 2 = 1 




5 10 15 

x, lattice units 



I 



to 




5 10 15 

x, lattice units 



FIG. 1. Contour plot of transverse anisotropy E of Eq. (g) 
created by: a - one edge dislocation along the y-axis; b - 
randomized array of dislocations. (Grey scales are arbitrary.) 



Due to the terms of the type S X S Z , etc., in the Hamil- 
tonian, the local easy axis deviates from the z-direction. 
In a locally rotated coordinate system (x 1 , y', z') that re- 
stores the normal form of the crystal field, saving terms 
linear on deformation, one obtains 



n = -DSl - H Z S Z 



E{Sl 



fla'Sx', (5) 



where 



E = giD(e xx - e y 



:{gz£xz + 9iU xz )H z . (6) 



A few observations are in order. Dislocation generates 
the transverse anisotropy of strength E that decays as 
l/r on the distance from the dislocation core. Due to 
the slow decay of E, a single dislocation produces spin 
tunneling at a large number of crystal sites. The field ap- 
plied along the c-axis of the crystal generates the trans- 
verse field due to the local rotation of the easy axis by the 
dislocation. Above the XY plane (at z > 0) the trans- 
verse field is directed along the hard axis, while at z < 
the transverse field is along the medium axis. 

For numerical work we choose gi = g% = .93 = .94 = 1 
(see, e.g., Refs. [jll| and references therein). The spatial 
dependence of the transverse anisotropy E due to a single 
edge dislocation, described by Eqs. (|J) and (||), is shown 




FIG. 2. Relaxation curves for the 100 x 100 Mni2 lattice 
around a single edge dislocation along the j/-axis at T — 2 K 
and H — 0. Faster relaxation of Mni2 molecules caused by the 
dislocation leads to a nonexponential form of R(t/r), where 
t is the relaxation time for an ideal sample. Dashed line: 
Pure thermal exponential in the ideal Mni2 crystal. Inset: 
Tunneling levels mi, = 0,-1, —2, —3, —4 at different points of 
the lattice (mt, = in the small lightmost regions corresponds 
to the 65K barrier in an ideal crystal). 



in Fig. [l]a. The effect of many dislocations is additive 
as long as the linear elastic theory is applied, which is 
correct outside dislocation cores. Fig. [lb shows a typi- 
cal pattern of the transverse anisotropy produced by an 
array of edge dislocations, obtained by randomization of 
the quadratic lattice of dislocations with alternating ori- 
entation of the extra plane. 

Due to the dislocations the local energy barrier be- 
tween spin-up and spin-down states is lower than in the 
ideal crystal. It can be calculated perturbatively from 
Eq. (||) . In addition to the barrier reduction for a classical 
spin, the barrier in the quantum case is further lowered 
at the discrete values of H z , 



H zk = k^D* - E 2 , fc = 0,±l,...,±2S (7) 
that provide resonant tunneling between matching spin 



levels m and m' = —m — k P,p[,p|j2ll 
effective height of the barrier, 



At H 7 = H zk the 
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is determined by m — mj < that corresponds to the 
lowest pair of levels, m and m', whose tunneling splitting, 
A mm /, is greater than the sum of their widths, T mm i = 
r m + r m / . Eq. (||) , even though it accounts for only one 
tunneling pair of levels, is a good approximation for the 
relaxation problem (see discussion after Eq. (5.15) and 
Fig. 5 of Ref. ||). The splitting is determined by high 
powers of E and H x >, and, thus, strongly depends on 
coordinates [see Eq. (2.6) of Ref. [|| for the effect of the 
transverse field and Eq. (4) of Ref. j22) for the effect of 
transverse anisotropy]. Consequently, the dependence of 



m,b on T mm i is weak. In the kelvin range we will use the 
experimental value T mm > ~ 200 Oe |p3[ . 

To illustrate the barrier reduction by dislocations in 
the thermally activated regime (kelvin range), we nu- 
merically diagonalized Eq. (||) for the 100 x 100 Mni2 
lattice around a single edge dislocation at T = 2 K and 
H = 0. Using the criterion described above, we deter- 
mine tunneling levels mf, and the effective barriers C/? ff 
from Eq. (||) for all sites i of the pattern and construct 
the relaxation curve 



R(t) 



N ^ 



exp(-a^/r), 



(9) 



where a* 



DS 2 is the unperturbed 



barrier, and r = TQe u °°/ T is the relaxation time far from 
the dislocation. One can see from Fig. || that the effect 
of dislocations on relaxation is profound. Even for one 
dislocation per 100 x 100 sites in the XZ plane, which cor- 
responds to the concentration of dislocations as small as 
c = 10~ 4 , more than half of Mni2 molecules relax faster 
than in the ideal crystal. The list of energy barriers reads 
65K (64), 64.35K (4700), 62.38K (4944), 59.06K (244), 
54.31K (36), 48.32K (8), 30.73K (4), with the number of 
corresponding crystal sites shown in the brackets. 

At temperatures below IK, where spin relaxation oc- 
curs via tunneling from the ground state, the influence 
of dislocations becomes even more dramatic. In contrast 
to the small barrier reduction in the thermally activated 
regime, here the role of dislocations is to provide the 
main source of spin relaxation. The resulting level split- 
tings are distributed over many decades. Fig. ||a shows 
distributions of the ground-state level splittings obtained 
numerically for a randomized array of edge dislocations 
with concentration c = 10~ 2 for k = and k — 8 reso- 
nances. Since outside dislocation cores the perturbations 
of the uniaxial Hamiltonian are small, one can scale the 
splitting distributions for even values of k using the per- 
turbative Eq. (4) of Ref. |22| for the ground-state split- 
ting (m = —5, ml = —m — k = S — k), 
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Then the distribution of ln(A ki /gk)/ (S-k/2) = Ei/(8D) 
does not depend on fc, see Fig. For k =/= there is 
also a transverse field H x i in addition to the transverse 
anisotropy E in Eq. (||), which should modify Eq. (10). 
This modification, however, is small away from the dislo- 
cation cores. A more significant effect of the transverse 
field generated by the dislocations is unfreezing of tun- 
neling resonances with odd values of k. Our calculation 
shows that the distribution of the level splittings for odd 
k is similar to that for even k in Fig. ||. 

Let us study now spin relaxation produced by sweeping 
the field H z through a tunneling resonance. The popu- 
lation of the metastable state with m — —S obeys [p| 
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FIG. 3. a - Distribution of level splittings (in kelvin) cre- 
ated by a randomized array of edge dislocations with c = 10~ 2 
for resonances with k — and k — 8; b - scaled representation 
with g k of Eq. @. 
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where s~s ~ £s-fc = Vkt and Vk = (2S — k)dH z /dt. Inte- 
grating this equation one obtains the probability for the 
molecule to remain in the metastable well after crossing 
the resonance, P = exp(— 7rA 2 /ufc), which resembles the 
Landau-Zener formula (see also Rcf. [0). The relaxation 
curve R(vk) for a sample with dislocations 



R{v k ) 



iv^> xp( - 
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at different resonances fc is shown in Fig. At odd k 
relaxation is produced by the transverse field H x i in com- 
bination with the transverse anisotropy. For k — 1 the 
effect of H x > is still small and the corresponding relax- 
ation curve is noticeably shifted to the right. 

Due to the wide distribution of splittings, the sweep- 
ing rate needed to make different Mni2 molecules to re- 
lax stretches over many decades. Consequently, R(vk) 
can only be plotted on the log scale. On that scale 
any individual exponential becomes a step function, 



1.0- 



0.8- 



0.6- 



0.4- 





i , i 






\ \ \2\ 
\ \ \ 3 W 






\ \ 4\ \ 






\6\ \\ 




a \ \ 


\ *\ \ \ 




k = 10 \ 






Experimental window , 
-1 1 1 1 1 1 1 r 


- 10^^^^ 
1 1 1 1 





10 15 20 25 30 




8 10 

ln[g t (7i/^)" 2 ]/(5-^/2) 
FIG. 4. a - Spin relaxation as a function of the field sweep 
rate Vk = (25* — k)dH z /dt for a randomized array of edge 
dislocations with c = 10~ 2 at different tunneling resonances 
k. b - scaled representation with of Eq. (0). dH z /dt is 
measured in K 2 ~ 10 11 T/s. 



exp(— 7rAf /vk) => 0(1 — 7rAf /ufc). Since i?(ufc) represents 
the fraction of Mni2 molecules which have 7rA^ < Vk 
and, thus, do not relax, R(vk) and the distribution of 
splittings are related as the integral and the derivative. 
Correspondingly, R(vk) for different k can be scaled the 
same way as in Fig. [3|b. Deviations from the perfect scal- 
ing in Fig. [|b are due to the transverse field. 

The effect of dislocations, even at moderate concen- 
trations, appears to be much stronger than the effect of 
transverse fields Hj_ from nuclear spins and dipole in- 
teractions. Simple arguments of the perturbation theory 
show that the effect of dislocations becomes comparable 
with that of H x at E ~ H\/{DS 2 ). With the help of 
Eqs. (|) and @ one finds that for c > 1CT 6 the effect 
of dislocations on tunneling is greater than the effect of 
hyperfine and dipole fields at almost all crystal sites. 

We have shown that dislocations in Mni2 crystals 
should be the main source of spin tunneling in the kelvin 
and subkelvin temperature range. Local rotations of the 
easy axes due to dislocations give rise to the effective 
transverse field which unfreezes odd resonances. Tunnel- 



ing lifetimes due to dislocations are spread over many 
decades, which results in the stretched relaxation, espe- 
cially pronounced in the subkelvin range. Distribution of 
the tunneling rates and the relaxation law obey scaling 
which does not depend on the type of crystal defect. That 
scaling should be seen in experiment. Eventually, it may 
also become possible to observe magnetization patterns 
on the surface of Mni2 crystals, shown in Fig. [TJ. 

We thank Jonathan Friedman for useful remarks. This 
work has been supported by the NSF Grant No. 9978882. 
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